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Optical turbulence is the term that is used being based on the analogy between hydrodynamic 
and optical equations. The phenomenon of turbulent photon filamentation occurs in lasers and 
other active optical media at high Fresnel numbers. A description of this phenomenon is suggested. 
The solutions to evolution equations are presented in the form of a bunch of filaments chaotically 
distributed in space and having different radii. The probability distribution of patterns is defined 
characterizing the probabilistic weight of different filaments. The most probable filament radius and 
filament number are found, being in good agreement with experiment. 
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1. Introduction 



Different physical systems often display interesting similarity in their behaviour. This concerns, for example, 
hydrodynamic and optical systems. There exists a well established analogy between the evolution equations describing 
optical phenomena and hydrodynamic equations for compressible viscous liquid [1-3]. The Presnel number for optical 
systems plays the same role as the Reynolds number for fluids. In optical systems, with increasing Fresnel number, 
there appears chaotic behaviour similar to fluid turbulence arising with increasing Reynolds number. Because of 
this analogy, one terms the optical turbulence [1-3] the chaotic phenomena occurring in high-Fresnel-number optical 
systems. 

For an optical cavity of the standard cylindrical shape of radius R and length L, the Fresnel number is F = ttR^/XL, 
with A being the optical wavelength. Physical processes accompanying the route from the regular behaviour to the 
turbulent regime, arising with increasing Fresnel numbers, are practically the same in diflferent types of active optical 
media, such as lasers and photorefractive crystals. 

When Fresnel numbers are very small, _F <C 1, there exists the sole transverse mode almost uniformly flUing the 
optical cavity. At Fresnel numbers around F ~ 1, several transverse modes emerge seen as bright spots in the 
transverse cross-section. These spatial structures are regular both in space and in time, since in space they form 
ordered geometric arrays, as polygons, and because in time they are either stationary or oscillating periodically. 
These transverse structures are imposed by the cavity geometry and correspond to the empty-cavity Gauss-Laguerre 
modes. The appearance of such regular geometric structures is well understood theoretically, with the theory being 
in good agreement with experiments both for lasers [4-9] as well as for photorefractive crystals [10-12]. For Fresnel 
numbers up to F « 5, the number of regular modes is proportional to F^. 

The transition from the regular behaviour to a turbulent state occurs around the Fresnel number F « 10. For 
Fresnel numbers F > 15, the arising spatial structures are principally different from those defined by the empty-cavity 
Gauss-Laguerre modes. The related modal description becomes no longer relevant and the boundary conditions do 
not play the role any more. In the medium, there appears a large number of filaments stretched along the cavity 
axis. In the transverse cross-section, they are distributed chaotically. The correlation length is much shorter than the 
mean distance between filaments, so that they radiate independently of each other, aperiodically flashing in time. The 
number of filaments, contrary to the regular regime, is proportional to F. This chaotic filamentation has been observed 
in both photorefractive crystals [10-12], e.g. in Bii2Si02o, as well as in many lasers, such as Ne, Tl, Pb, N2, and 
vapor lasers [13-17], CO2 and dye lasers [18-26]. The specific random properties of the observed filamentation, such 
as chaotic distribution of filaments in space, the absence of correlations between them, and their aperiodic fiashing, 
are analogous to the properties of turbulent fluids [1-3]. Because of this analogy of the chaotic photon filamentation in 
high-Fresnel-number optical media with turbulent phenomena in high-Reynolds-number viscous fluids, one commonly 
calls [18-26] the former the optical turbulence or, because of the formation of bright filaments with a high density of 
photons, it can be called the turbulent photon filamentation. 

Contrary to the low-Fresnel number case, with regular optical structures, for which a good theoretical understanding 
has been achieved [4-12], for the high-Fresnel-number regime of the turbulent photon filamentation, there has been 
no persuasive theory suggested. This problem was addressed in Refs. [27-30]. However, because of the complexity of 
the phenomenon, only some oversimplified stationary models were considered. In the present paper, the description 
of the phenomenon of turbulent photon filamentation is suggested, based on realistic evolution equations. It has 
become possible to solve now this problem owing to the recently developed method for treating nonlinear differential 
equations, called the scale separation approach [31-33]. 

2. Turbulent photon filamentation 

For the system of N resonant atoms interacting with electromagnetic field, we use the standard dipole approxima- 
tion. The evolution equations to be considered are those for the averages of atomic operators, 



where a~ is a transition operator and is the population-difference operator. Equations for the quantities (1) can be 

easily obtained following the usual way, by writing the Heisenberg equations of motion, eliminating the field variables, 
and employing the standard semiclassical and Born approximations. All this machinery is well known and can be 
found, for instance, in the book [34]. In order to present resulting equations in a compact form, it is convenient 
to introduce the notation f = fo + frad for an effective field acting on an atom. This field consists of the term 
/o = —id ■ Fo, due to the cavity seed field, and of the term 



u{f, t) = < a {f,t) > , 



s{f,t) = < a''{r,t) > 



(1) 




(2) 
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describing the effective interaction of atoms through the common radiation field. Here p is the density of atoms, 
d = doCd is the transition dipole, and ip{r) = exp(ifco|r|)/fcoH) with ko = coq/c and 7 = 4A:Qdo/3. The cavity seed 
field selecting a longitudinal propagating mode, is taken, as usual, in the form of plane waves. 



The resulting equations are 



_= _(ia;o+72)n + 5/, -= -2(«*/ + r«)-7i(s-C) 



d\ 



2 



dt 



= -2j2\u\' + s{u*f + ru), (3) 



where wq is the transition frequency, 71 and 72 are the longitudinal and transverse relaxation widths, respectively, 
and C is the pumping parameter. 

We consider a cylindrical cavity, typical of lasers, with radius R and length L. For the case of high Fresnel numbers 
F 10, the cavity can house several transverse modes. Therefore, it is reasonable to look for solutions of Eqs. (3) 
in the form of a bunch of Nf filaments, so that 

Nf Nf 

u{r,t) = Un{r,t)Qn{x,y) e''=^ , s{r,t) = ^ s„(f, i)e„(a;, y) , (4) 

n=l n=l 



where 



e„(x, y) = Q{hl -{x- x,,f + (y - y^f) 



is a unit-step function; the pair Xn and yn defines the axis of a filament: and 6„ is the radial distance from the 
filament axis, at which the solutions decrease by an order of magnitude. I f the f ilament profile is approximately of 
normal law exp(— r^/2r^), with r„ being the filament radius, then &„ = \/2 In 10 r„. The filaments do not interact 
with each other, because of which they cannot form a regular space structure, and the locations of their axes in the 
transverse cross-section are random. The absence of interactions between filaments is due to the fact that the kernel 
(^(r), playing the role of a Green function in the effective interaction field (2), fastly oscillates and diminishes with 
increasing r. Let us stress that the presentation of solutions in the form of a bunch of filaments (4) does not exclude 
the possible case of a sole filament almost uniformly filling the whole sample. This is because the number of filaments 
and their radii will be defined in a self-consistent way. Substituting the form (4) into Eqs. (3), we obtain a system of 
equations for w„, s„, and \un\^ ■ To simplify these equations, we employ the mean-field approximation for the averages 

U{t)= / Un{f,t)df, S{t)= / Sn{f,t)df, (5) 

and for the similarly defined where the integration is over the volume Vn = T^b^L of a cylinder enveloping a 

filament. For the simplicity of notation, the index n labelling filaments is dropped from the left-hand sides of Eqs. 
(5). In this way, Eqs. (3) are transformed to the system of equations 

fill -* 

— = + T)u - isd ■ El e-"^* , 

J = -4gj2\uf - 71 (s - C) - 4 Im (u*d- E^ e"-*) , (6) 



^ = -2r|w|2 + 2s Im (u*d- Er e"''^*) 



for the functions (5). Here the notations for the collective frequency Q = + g'l2S and for the collective width 
r = 72(1 — gs) are used, with the effective coupling parameters 

g= 37P j sin[fco|f - fj -k{z - z')] 



472V„ Jv^ ko\f-f'\ 
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a'= f cos[ko\f-r'\-kiz-z')] 

Equations (6) describe the evolution of a filament. These evolution equations can be analysed by using the scale 
separation approach [31-33]. For applying the latter, we, first, take into account the existence of usual small parameters 
71 /wo ^ 1 J^nd 72/i^'o 1, and also wc consider the quasircsonancc case, when |A|/ciJo ^ with the detuning 
A = w — (jjQ. Then the functional variable u in Eqs. (6) is to be classified as fast, as compared to the slow variables s 
and The latter are quasi-invariants for the fast variable. Following the scale separation approach [31-33], we solve 
the equation for the fast variable u, which is straightforward when s is a quasi-invariant. The found solution for u 
has to be substituted into the right-hand sides of the equations for the slow variables s and j-uj^, with averaging these 
right-hand sides over explicitly entering time of fast oscillations. This procedure, complimented by the transformation 

I |2 2 M ■ -^1 P /^\ 

where a is assumed to be small, a <^1, leads to the equations 

— = -4572M; -71 (s-g , —= -272(1 - ffsjw (9) 

for the slow variables. 

The properties of solutions to Eqs. (9) essentially depend on the value of the effective coupling g defined in Eq. 
(7). From this definition, it is evident that the value of g depends on the characteristics of the considered filament, 
in particular, on the filament radius r„ which enters Eq. (7) through the radius 6„ of the enveloping cylinder. It 
is convenient to introduce the dimensionless parameter (3 = irb^/XL and to consider g — g{(3) as a function of this 
parameter in the domain < P < F. Thus, each filament can be characterized by its radius or, equivalently, by the 
related parameter /?. The bunch of filaments parametrized by (3 forms the solutions (4). The stability properties of 
filaments with different values of f3 are different, which can be described by defining the probability distribution of 
patterns corresponding to filaments with (3 varying in the interval [0, -F]. The latter can be done as follows. Any 
system of evolution equations, like Eqs. (9), can always be presented in the normal form dy/dt = v, in which y = {yi}, 
and V = {vi} is a, velocity field. In the case of Eqs. (9), we have yi = s, j/2 = w. 

As is known from statistical mechanics, a probability p is connected with entropy S by the relation p ^ . In 
the nonequilibrium case, it is more appropriate to count the entropy S{t), which is a function of time, from its initial 
value S{0), thus, considering the entropy variation AS{t) = S{t) — S{0). The entropy is defined as the logarithm of 
an elementary phase volume, S{t) = In |(5r(i)|, where, in the nonequilibrium case, ST{t) = Ylj 5yj{t). The elementary 
phase volume can be expressed as 5T{t) — Hi -^ij(^)^2/i(0) through the elements Mij{t) = Syi{t)/Syj{0) of the 
multiplier matrix M{t) = [Mij{t)]. Hence, the entropy variation is AS{t) = TrL, with the matrix L = [Lij] being 
composed of the elements L^j = In |My |. By varying the evolution equations, one gets the equation dM/dt = JM for 
the multiplier matrix M{t), where J = [Jij] is the Jacobian matrix with the elements Jij = Svi/Syj. In this way, the 
probability p ~ e"'^'^ acquires the form p ~ exp{— TrL}. From the equation for the multiplier matrix, one has the 
entropy variation AS{t) = /g K{t')dt', where K = TrJ called the contraction rate [35]. Therefore, for the probability 
distribution of patterns, labelled by a parameter /?, one obtains 







piP,t) = exp<i - / K{p,t') dt' \ , (10) 



where 

Z{t)= I exp<;- / K{f3,t')dt') dp (11) 



t 







is the normalization factor. As is clear, that pattern is preferred over the other which possesses a higher probability 
weight. This implies, because of the form (10), that the most preferable pattern is that for which the local contraction 

1 /■* 

K{t) = - / K{t')dt' (12) 
* Jo 

is minimal. At the initial stage of pattern formation, one may write p ~ exp{— i(r(/3, 0)t}. Hence, the most probable 
pattern to be formed is that corresponding to the minimal contraction rate K{(3, 0). 
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Considering the evolution equations (9), we keep in mind the standard way of exciting laser media by means of a 
nonresonant pumping described by the pumping parameter C > 0. No resonant fields are involved, so that s(0) < 0. 
It is easy to calculate the contraction rate K{f3, 0) corresponding to Eqs. (9), which is 

0) = -71- 272(1 -5So), (13) 

where sq = s(0). The minimum of this contraction rate is provided by the maximum of <? = g{P). The maximum of 
g{l3), with g{l3) given by Eq. (7), occurs at /? = 0.96. Prom here we find the most probable filament radius 

rf = 0.26VAL . (14) 
The most probable number of filaments can be evaluated from the normalization integral 

J s{f,t)df= C , 

where the integration runs over the whole volume of the sample, V = ttR^L. Considering this normalization integral 
for the stage when the filaments have already been formed and the population difference inside a filament of radius 
r/ has reached a value close to (, we obtain 

^/=(^)'=4.71F. (15) 

The number of filaments is proportional to the Presnel number F, which is common for the turbulent photon filamen- 
tation. 

Calculating the filament radius r/ for the variety of experiments with different lasers, we find that the theoretical 
values of r/ are in perfect agreement with all experimental data available. Thus, for Ne, Tl, Pb, N2, and laser 
[13-17], we find r/ « 0.01 cm and for CO2 and dye lasers [18-24], we have r/ « 0.08 cm and r/ « 0.01 cm, respectively. 
It is interesting that the coefficient 4.71 in the dependence (15) for the number of filaments is also in good agreement 
with those experiments where it was measured [10-12]. 

3. Conclusion 

The solution for the problem of turbulent photon filamentation occurring in resonant media at high Fresnel numbers 
is suggested. The consideration is based on realistic evolution equations for resonant atoms under conditions typical 
of lasers. The solutions to these equations have the form corresponding to a bunch of filaments with different radii. 
The probability distribution of filament radii is found self-consistently from the evolution equations. The results are 
in good agreement with all experiments on the turbulent photon filamentation in laser media [13-24]. The method 
presented in this paper can also be employed for considering other physical systems with arising spatio-temporal 
structures, when one confronts the so-called problem of pattern selection. 

This work has been supported by the Bogolubov-Infeld Grant of the State Agency for Atomic Energy, Poland and 
by a Grant of the Sao Paulo State Research Foundation, Brazil. 
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